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Abstract 



The functional significance of correlations between action potentials of neurons is still a matter of vivid 
debates. In particular it is presently unclear how much synchrony is caused by afferent synchronized events 
and how much is intrinsic due to the connectivity structure of cortex. The available analytical approaches 
based on the diffusion approximation do not allow to model spike synchrony, preventing a thorough analysis. 
Here we theoretically investigate to what extent common synaptic afferents and synchronized inputs each 
contribute to closely time-locked spiking activity of pairs of neurons. We employ direct simulation and extend 
earlier analytical methods based on the diffusion approximation to pulse-coupling, allowing us to introduce 
precisely timed correlations in the spiking activity of the synaptic afferents. We investigate the transmission 
of correlated synaptic input currents by pairs of integrate-and-fire model neurons, so that the same input 
covariance can be realized by common inputs or by spiking synchrony. We identify two distinct regimes: In the 
limit of low correlation linear perturbation theory accurately determines the correlation transmission coefficient, 
which is typically smaller than unity, but increases sensitively even for weakly synchronous inputs. In the limit 
of high afferent correlation, in the presence of synchrony a qualitatively new picture arises. As the non-linear 
neuronal response becomes dominant, the output correlation becomes higher than the total correlation in the 
input. This transmission coefficient larger unity is a direct consequence of non-linear neural processing in the 
presence of noise, elucidating how synchrony-coded signals benefit from these generic properties present in 
cortical networks. 

Author summary 

Whether spike timing conveys information in cortical networks or whether the firing rate alone is sufficient 
is a matter of controversial debates, touching the fundamental question how the brain processes, stores, and 
conveys information. If the firing rate alone is the decisive signal used in the brain, correlations between action 
potentials are just an epiphenomenon of high convergence and divergence the cortex' connectivity, where pairs 
of neurons share a considerable fraction of common afferents. Due to the membrane leakage, the small synaptic 
amplitudes and the non-linear threshold, nerve cells exhibit lossy transmission of correlation, if the correlation 
originates from shared synaptic inputs. However, the membrane potential of cortical neurons often displays 
non-Gaussian fluctuations, a hallmark of synchronized synaptic afferents. Moreover, synchronously active 
neurons have been found to reflect behavior in primates. In this work we therefore contrast the transmission of 
correlation due to shared afferents and due to synchronously arriving synaptic impulses for leaky neuron models. 
We not only find that neurons are highly sensitive to synchronous afferents, but that they are generically able 
to perform noise suppression on synchrony coded signals, a computational advantage over rate signals. 

1 Introduction 

Simultaneously recording the activity of multiple neurons provides a unique tool to observe the activity in the 
brain. The immediately arising question of the meaning of the observed correlated activity between different 
cells [1, 2] is tightly linked to the problem how information is represented and processed by the brain. This 
problem is matter of an ongoing debate [3] and has lead to two opposing views. In one view, the high 
variability of the neuronal response [4] to presented stimuli and the sensitivity of network activity to the 
exact timing of spikes [5] suggests that the slowly varying rate of action potentials carries the information 
in the cortex. A downstream neuron can read out the information by pooling a sufficient number of merely 
independent stochastic source signals. Correlations between neurons may either decrease the signal-to-noise 
ratio of population signals [6] or they may enhance the information in such population, depending on the 
readout mechanism [7]. Correlations are an unavoidable consequence of cortical connectivity where pairs of 
neurons share a considerable amount of common synaptic afferents [8]. Recent works have reported very 
low average correlations in cortical networks on long time scales [9], explainable by an active mechanism of 
decorrelation [10, 11, 12]. On top of these correlations inherent to cortex due to its connectivity, a common 
and slowly varying stimulus can evoke correlations on a long time scale. 

In the other view, on the contrary, theoretical considerations [13, 14, 15, 16] argue for the benefit of 
precisely timed action potentials to convey and process information by binding elementary representations into 
larger percepts. Indeed, in frontal cortex of macaque, correlated firing has been observed to be modulated in 
response to behavioral events, independent of the neurons' firing rate [17]. On a fine temporal scale, synchrony 
of action potentials [18, 19, 20] has been found to dynamically change in time in relation to behavior in primary 
visual cortex [21] and in motor cortex [17, 22]. The observation that nearby neurons exclusively show positive 
correlations suggests common synaptic afferents to be involved in the modulation of correlations [23]. In this 
view, the measure of interest are correlations on a short temporal scale, often referred to as synchrony. 

The role of correlations entails the question whether cortical neurons operate as integrators or as coincidence 
detectors [18, 24]. Recent studies have shown that single neurons may operate in both regimes [25]. If the 
firing rate is the decisive signal, integrator properties become important, as neural firing is driven by the 
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mean input. As activity is modulated by the slowly varying signal, correlations extend to long time scales 
due to co-modulation of the rate. Integrators are thus tailored to the processing of rate coded signals and 
they transmit temporal patterns only unreliably. Coincidence detectors preferentially fire due to synchronously 
arriving input. The subthreshold membrane potential fluctuations reflect the statistics of the summed synaptic 
input [26], which can be used to identify temporally precise repetition of network activity [27]. A direct probe 
for the existence of synchronous activity are the resulting strong deflections due to synchronous arrival of 
synaptic impulses. Such non-Gaussian fluctuations have indeed been observed in auditory cortex in vivo [28] 
and in the barrel cortex of behaving mice [29]. In this regime, coincidence detector properties become crucial. 
Coincidence detectors are additionally sensitive to stimulus variance [25, 30] and exhibit correlations arising 
from precisely timed firing. This type of correlation is unaffected by firing rate, can encode stimulus properties 
independently and moreover arises on short time scales [25]. 

The pivotal role of correlations distinguishing the two opposing views suggests to ask the following question: 
Can the experimentally observed synchrony between the activity of two neurons be explained by common input 
or is synchrony in the input required? If common input is sufficient, synchrony is just a side effect of rate 
coding. However, if synchrony in the input is required, this synchrony is likely to propagate information 
though the network, as it appears at task-specific times [17, 21, 22]. A functional interpretation is assigned to 
synchrony by the picture of the cell assembly [13, 14, 31, 32], where jointly firing neurons dynamically form 
a functionally relevant subnetwork. Due to the local connectivity with high divergence and convergence, any 
pair of neurons shares a certain amount of input. This common input may furthermore exhibit spike synchrony, 
representing the coherent activity of the other members of the cell assembly. In the assembly picture, the 
synchronous input from peer neurons of the same assembly is thus considered conveying the signal, while the 
input from neurons outside of the assembly is considered as noise [33]. 

One particular measure for assessing the transmission of correlation by a pair of the neurons is the trans- 
mission coefficient, i.e. the ratio of output to input correlation. When studying spiking neuron models, the 
synaptic input is typically modeled as Gaussian white noise, e.g. by applying the diffusion approximation to 
the leaky integrate-and-fire model [34]. In the diffusion limit, the transmission coefficient of a pair of model 
neurons receiving correlated input mainly depends on the firing rate of the neurons alone [35, 36]. For low 
correlations, linear perturbation theory well describes the transmission coefficient, which is always below unity, 
i.e. the output correlation is bounded by the input correlation, pairs of neurons always lose correlation. Ana- 
lytically tractable approximations of leaky integrate-and-fire neural dynamics have related the low correlation 
transmission to the limited memory of the membrane voltage [37]. The transmission is lowest if neurons are 
driven by excitation and inhibition, when fluctuations dominate the firing. In the mean driven regime the 
transmission coefficient can reach unity for integral measures of correlation [37], 

Understanding the influence of synchrony among the inputs on the correlation transmission requires to 
extend the above mentioned methods, as Gaussian fluctuating input does not allow to represent individual 
synaptic events, not to mention synchrony. Therefore, in this work we introduce an input model that extends 
the commonly investigated Gaussian white noise model. We employ the multiple interaction process (MIP) 
[38] to generate an input ensemble of Poisson spike trains with a predefined pairwise correlation coefficient. 
We use these processes containing spike synchrony as the input common to both neurons and model the 
remaining afferents as independent Poisson spike trains. Furthermore, contrary to studies that measure the 
integrated output correlation (count correlation) [35, 36], we primarily consider the output correlation on the 
time scale of milliseconds, i.e. the type of correlation determined by the coincidence detection properties of 
neurons. 

In Sec. 2 we first introduce the neuron and input models. In Sec. 2.1 we study the impact of input 
synchrony on the firing properties of a pair of leaky integrate-and-fire neurons with current based synapses. 
Isolating and controlling this impact allows us to separately study the effect of input synchrony on the one hand 
and common input on the other hand on the correlation transmission. In Sec. 2.3 and Sec. 2.4 we present a 
quantitative explanation of the mechanisms involved in correlation transmission, in the limit of low and high 
correlation, respectively, and show how the transmission coefficient can exceed unity in the latter case. In 
Sec. 3 we summarize our findings in the light of previous research, provide a simplified model that enables an 
intuitive understanding and illustrates the generality of our findings. Finally, we discuss the limitations of our 
theory and consider possible further directions. 

2 Results 

The neuronal dynamics considered in this work follows the leaky integrate-and-fire model, whose membrane 
potential V(t) obeys the differential equation 

V(t) <- V r if V(t) > Vo, 
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where r m is the membrane time constant, Vq the resting potential, V$ the firing threshold, and V r the reset 
potential of the neuron. The neuron is driven by excitatory and inhibitory afferent spike trains s cxc (t) — 
wJ2j$(t — *cxc) anc ' s inh(£) = — 9 ' w J2k^(^ ~ ^Inh) where w ls the excitatory synaptic weight and P exc 
and tf nh are the arrival time points of excitatory and inhibitory synaptic events, respectively. s oxc /i nn denote 
the weighted sum of all afferent excitatory and inhibitory events, respectively. Inhibitory events are further 
weighted by the factor —g. Each single incoming excitatory or inhibitory event causes a jump of the membrane 
potential by the synaptic weight w or -gw, respectively, according to (1). Whenever the membrane potential 
reaches the threshold V$ the neuron fires a spike and the membrane potential is reset to V r after which it is 
clamped to that voltage for a refractory period of duration t t . In the current work we measure the correlation 
between two spike trains si and Sj on the time scale t as 

^ / CovK.n?] \ 
\^Var[nJ-]Var[nJ]/ r 

where nj is the spike count of spike train Si in a time window r and the average ()t is performed over the 
T/t time bins of a stationary trial. Except for Sec. 2.3 where we choose r = 100 ms, throughout our work 
we choose r = 1 ms and refer to Po U " ls as output spike synchrony. 

We investigate the correlation transmission of a pair of neurons considering the following input scenario. 
Each neuron receives input from N presynaptic neurons of which fN are excitatory and (l — f)N are inhibitory. 
Both neurons share a fraction c € [0, 1] of their excitatory and inhibitory afferents. Fig. 1A shows a schematic 
representation of the input to neurons i = 1,2. Each source individually obeys Poisson statistics with rate v- m . 
We assume that the two neurons under consideration are part of a subnetwork performing some function. They 
receive a number cfN common excitatory inputs from neurons belonging to the same functional unit. We 
want to investigate how synchrony can be used to convey information within the functional subnetwork. We 
therefore allow synchronization of spikes to occur in the cfN common excitatory inputs, whereas the (l — c)fN 
excitatory disjoint afferents are assumed to originate from independent source networks. The assumption of 
independence also holds for all (1 — f)N inhibitory sources. Here in particular, we use a multiple interaction 
process (MIP) [38] to model synchronous spike events among the common, excitatory afferents. In this model 
each event of a mother Poisson process of rate A m is copied independently to any of the cfN child spike 
trains with probability p, resulting in a pairwise correlation coefficient of p between two child spike trains. 
By choosing the rate of the mother spike train as X m = ^ the rate of a single child spike train is v- m and 
independent of p. 

Fig. IB, C, D and E show that the amount of pairwise correlations in the common input has a strong 
impact on the variance and correlation of the free membrane potentials {a 2 ,p ln ) and therefore on the output 
firing rate and output spike synchrony {v OVLt , p\™). Let us first consider the case of p= 0, i.e. the absence 
of synchronous events in the input. As expected, the free membrane potential variance a 2 remains constant 
throughout the whole range of c, as does the firing rate v out (Fig. IB and D). Figure Fig. 1C shows the 
correlation of the free membrane potential of a neuron pair, normalized by the common input fraction c. 
As expected, for p = the input correlation is only determined by the common input fraction and thus 
Pin = c. Hence, the output synchrony observed for p = in Fig. IE is solely due to the correlation caused 
by common input and describes the often reported correlation transmission function of the integrate-and-fire 
model [35, 36], where for < c < 1 the output spike synchrony is always well below the identity line, which 
is in full agreement with the work of [35]. 

Let us now consider the case of p > 0. In Fig. IB and D we observe that even small amounts of input 
synchrony result in an increased variance of the free membrane potential, which is accompanied by an increase 
of the output firing rate. While for weak input synchrony the increase of a and v ont is only moderate, in 
the extreme case of strong input synchrony (p = 0.1) a becomes almost ten-fold higher and v out increases 
more than three-fold compared to the case of p = 0. Fig. 1C shows that input synchrony also has a strong 
impact on the correlation between the free membrane potentials of a neuron pair. For any p > the input 
correlation is most pronounced for high p and in the lower regime of c. Simulation results shown in Fig. IE 
suggest that this increase of input correlation is accompanied by an increased synchrony between the output 
spikes for p = 0.001 and p = 0.01. For strong input synchrony of p = 0.1 the output synchrony is always 
higher than the input correlation caused solely by the common input, except near c = and at c = 1. 

The output firing rates and output spike synchrony shown in figures Fig. ID and E bear a remarkable 
resemblance, most notably for lower values of c. Particularly salient is the course of these quantities for 
p = 0.1, which is almost identical over the whole range of c. These observations clearly corroborate findings 
from previous studies that have shown an increase of the correlation transmission of a pair of neurons with the 
firing rate of the neurons [35, 36]. Thus, we must presume that a substantial amount of the output synchrony 
observed in Fig. IE can be accounted for by the firing rate increase observed in Fig. ID. Furthermore, as 
Fig. 1C suggests, for any p > common input and the synchronous events both contribute to the correlation 
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Figure 1: A pair of integrate-and-fire model neurons driven by partially shared and correlated presynaptic 
events. A Each of the neurons i and j receives input from N sources, of which fN are excitatory and 
(1 — f)N are inhibitory. Both neurons share a fraction c of their excitatory and inhibitory sources, whereas 
the fraction (1 — c) is independent for each neuron. Schematically represented spike trains on the left of 
the diagram show the excitatory part of the input, the inhibitory input is only indicated. A single source 
emits spike events with a firing rate v ln , with marginal Poisson statistics. Correlated spiking is introduced 
in the cfN common excitatory sources to both neurons. This pairwise correlation is realized by means of a 
multiple interaction process (MIP) [38] that yields a correlation coefficient of p between any pairs of sources. 
In absence of a threshold, the summed input drives the membrane potential to a particular working point 
described by its mean fi and standard deviation a and the correlation coefficient p m = Cov[Vj, Vj]/(<Ji<jj) 
between the free membrane potentials Vi, Vj of both neurons. In presence of a threshold mean and variance 
of the membrane potential determine the output firing rate v out and their correlation in addition determines 
the output correlation p ut, calculated by (2). B-E Direct simulation was performed using different values of 
common input fraction c and four fixed values of input spike synchrony p (as denoted in E). Each combination 
of c and p was simulated for 100 seconds, gray coded data points show the average over 50 independent 
realizations. Remaining parameters are given in Table 1. Solid lines in B and C are calculated as (5) and (6), 
respectively. In C, for convenience, p[ n is normalized by the common input fraction c, so that p- m /c = 1 in 
absence of synchrony (p = 0). In E, we measure output spike synchrony p^ s - 
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Table 1: Parameters of the input and LIF neuron used in the simulations 



between the membrane potentials of a neuron pair. 

2.1 Understanding and isolating the effect of synchrony. These two observations - the increase of 
input correlation and output firing rate induced by input synchrony - foil our objective to understand the 
sole impact of synchronous input events on the correlation transmission of neurons. In the following we will 
therefore first provide a quantitative description of the effect of finite sized presynaptic events on the membrane 
potential dynamics and subsequently describe a way to isolate and control this effect. 

The synchronous arrival of k events has a fc-fold effect on the voltage due to the linear superposition 
of synaptic currents. The total synaptic input can hence be described by a sequence of time points V and 
independent and identically distributed (i.i.d) random number that assume a discrete set of synaptic 
amplitudes each with probability P(w' J ). The train of afferent impulses follows Poisson statistics with some 
rate A. Assuming small weights w and high, stationary input rate A, a Kramers-Moyal expansion [39, 40, 41] 
can be applied to (1) to obtain a Fokker-Planck equation for the membrane potential distribution p(V,i) 



Only the first two moments (if- 7 ) = Ylwi wJ P( w ^) ar| d ((w^) 2 ) = ~^2 w j (W) 2 P{w^) of the amplitude 
distribution enter the first (/i) and second (er 2 ) infinitesimal moments as [42, cf. Appendix Input-Output 
Correlation of an Integrate-and-Fire Neuron for a detailed derivation] 



H = \T m (wi)+V (4) 
v 2 = \r m \{{wif). 

In the absence of a threshold, the stationary density follows from the solution of S(V,t) = as a Gaussian 
with mean /i and variance a 2 . 

Equation (3) and (4) hold in general for excitatory events with i.i.d. random amplitudes arriving at Poisson 
time points. Given the K = cfN common excitatory afferents' activities are generated by a MIP process, the 
number of k synchronized afferents follows a binomial distribution P(k) = B(K,p,k) = (^)p fe (l — p) K ~ k , 
with moments (k) = Kp and (/c 2 ) = Kp(l — p + Kp). The total rate — (k) = v- ln K of arriving events is 
independent of p, as is the contribution to the mean membrane potential fj,. Further we assume the neurons 
to be contained in a network that is in the balanced state, i.e. g = f/(l — f), and that all afferents have the 
same rate u- m . Thus, excitation and inhibition cancel in the mean so that \i = Vq. Due to the independence 
of excitatory and inhibitory spike trains they contribute additively to the variance a 2 in (4). The variance due 
to (1 — f)N inhibitory afferents with rate v m is (1 — f)Nu- m g 2 F2, with F2 = \T m w 2 . An analog expression 
holds for the contribution of unsynchronized disjoint excitatory afferents. The contribution of K excitatory 
afferents from the MIP follows from (4) as £j» (k 2X j F2. So together we obtain 



a 2 = (cf(l-p + cfNp) + (l-c)f + g 2 (l- f)) Nv in F 2 (5) 
= (/(l ~cp + c 2 fNp) + g 2 (l - /)) Nv m Fi. 

Fig. IB shows that (5) is in good agreement with simulation results. We are further interested in describing 
the correlation p- m between the membrane potentials of both neurons. The covariance is caused by the 
contribution from shared excitation ^ (fc 2 ) F2, in addition to the contribution from shared inhibition c(l — 
f)Nvi n g 2 F2, which together result in the correlation coefficient 

p- m = (f(l-p + cfNp)+g 2 (l- f))cN Vin F 2 /a 2 . (6) 

Again, Fig. 1C shows that (6) is in good agreement with simulation results. In order to isolate and control 
the effect of the synchrony parameter p on (5) and (6), in the following we will contrast two distinct scenarios: 
In the first scenario, we generate a certain amount of input correlation p- m using common input alone without 



5 



A B 




■9 p i i i d u i i i l 

.25 .5 .75 1 .25 .5 .75 1 

input correlation Pin input correlation Pin 



Figure 2: Isolation and control of the effect of synchrony on the free membrane potential statistics. A,B 
Adjusted common input fraction c (A) and input firing rate v m (B) for different values of p (gray coded) 
that ensure the same variance and covariance as for p = 0. C Correlation coefficient p m normalized by c 
between the free membrane potential of a pair of neurons using the adjusted common input fraction c. D 
Standard deviation of the free membrane potential, using the adjusted firing rate v ln . The statistics of the 
free membrane potential measured in simulations in panels C and D are further verified via (6) and (5) (solid 
lines). 



any spiking synchrony p. In the second scenario we generate the same amount of input correlation p ln by 
using a certain amount of spike synchrony p > 0. In order to have the same working point in both cases, we 
need to keep the mean p and variance a 2 of the single neuron's membrane potential constant for a given fixed 
input correlation p m . At p = the input correlation p- m (6) equals the common input fraction c and does not 
depend on v- ln . For p > the same input correlation can be achieved by an appropriate decrease of c. To 
this end we solve (6) for c, obtaining the adjusted common input fraction c(pi n ,p) as the positive root of the 
resulting quadratic equation. 

Since c(pi n ,p) does not depend on the afferent rate we may adjust v m in order to keep the total variance 
(5) on a constant level. At p = and firing rate v\ n the variance is given by c 2 | p = (/ + g 2 (l — /)) Nv- m F2- 
For p > the same variance can be achieved by an appropriate decrease of v m . We therefore solve (5) for 
v m , obtaining the reduced firing rate Vi n (a 2 , p m: p). 

We evaluate this approach by simulating the free membrane potential of a pair of leaky integrate-and-fire 
neurons driven by correlated input. For different values of p, we choose c(p- ln ,p) and v- ln {a 2 , p lni p), shown 
in Fig. 2A and B, to keep the variance and the correlation constant. Fig. 2A shows that the adjustment of 
the common input fraction becomes substantial only for higher values of p: while for p = 0.001 the reduced 
c is only slightly smaller than c, for p = 0.1 and p m = 0.8 it is reduced to c = 0.21. Fig. 2B shows that 
even for small amounts of input synchrony, v- ln needs to be decreased considerably in order to prevent the 
increase of membrane potential variance (Fig. IB). In the extreme case of p- 1B = 1 and p = 0.1 (both neurons 
receive identical and strongly synchronous excitatory input) an initial input firing rate of 10 Hz needs to be 
decreased to v- ln = 0.15 Hz. Fig. 2C and D confirm that indeed both the correlation and the variance of the 
free membrane potential remain constant throughout the whole range of p ln and for all simulated values of p. 

2.2 Correlation Transmission. In order to study the transmission of correlation by a pair of neurons, we 
need to ensure that the single neuron's working point does not change with the correlation structure of the 
input. The diffusion approximation (3) suggests, that the decisive properties of the marginal input statistics 
are characterized by the first (p) and second moment (cr 2 ). As we supply balanced spiking activity to each 
neuron, the mean p is solely controlled by the resting potential Vq, as outlined above. For any given value 
of p and pi U , choosing the afferent rate Di n (a 2 , Pi n ,p) ensures a constant second moment a 2 . Consequently, 
Fig. 3B confirms that the fixed working point (p,a 2 ) results in an approximately constant neural firing rate 
^out for weak to moderate input synchrony p. For strong synchrony, fluctuations of the membrane potential 
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become non-Gaussian and the firing rate decreases; the diffusion approximation breaks down. 

In studies which investigate the effect of common input on the correlation transmission of neurons, the 
input correlation is identical to the common input fraction c [35, 36]. In the presence of input synchrony this 
is obviously not the case (Fig. 1C). Choosing the afferent rate and the common input fraction according to 
Pin (a 2 , pin, p) and c(pi n ,p), respectively, enables us to realize the same input correlation p ln with different 
contributions from shared inputs and synchronized events. We may thus investigate how the transmission of 
correlation by a neuron pair depends on the relative contribution of synchrony to the input correlation p- ln . 
Fig. 3A shows the output synchrony as a function of p ln for four fixed values of input synchrony p. As the 
input correlation is by construction the same for all values of p, changes in the output synchrony directly 
correspond to a different correlation transmission coefficient. Even weak spiking synchrony (p = 0.001) in 
the common input effectively increases the output synchrony, compared to the case where the same input 
correlation is exclusively caused by common input (p = 0). Stronger synchrony (p = 0.01 and p = 0.1) 
further increases this transmission. In Fig. 3B we confirm that the increase of output spike synchrony is not 
caused by an increase of the output firing rate of the neurons, but rather their rate remains constant up to 
intermediate values of p < 0.01. The drastic decrease of the output firing rate for p = 0.1 does not rebut 
our point, but rather strengthens it: correlation transmission is expected to decrease with lower firing rate 
[43, 36] for Gaussian inputs. However, here we observe the opposite effect in the case of strongly non-Gaussian 
inputs due to synchronous afferent spiking. We will discuss this issue in the following paragraph, deriving an 
analytical prediction for the correlation transmission. Moreover, we observe that the increased transmission is 
accompanied by a sharpening of the correlation function with respect to the case of p = (cf. Fig. 3C and 
D). 

For correlated inputs caused by common inputs alone (no synchrony, p = 0) or by weak spiking synchrony 
(p < 0.01) the transmission curves in Fig. 3A are always below the identity line. This means that the neural 
dynamics does not transmit the correlation perfectly, but rather causes a decorrelation. Recent work has shown 
that the finite life time of the memory stored in the membrane voltage of a leaky integrate-and-fire neuron 
is directly related to this decorrelation [37]. Quantitative approximations of this decorrelation by non-linear 
threshold units can be understood in the Gaussian white noise limit [35, 36]. For input correlation caused 
by spiking synchrony, however, we observe a qualitatively new feature here. In the presence of strong spiking 
synchrony (p = 0.1), in the regime of high input correlation (/?;„ > 0.8) the correlation transmission coefficient 
exceeds unity. In other words, the neurons correlate their spiking activity at a level that is higher than the 
correlation between their inputs. In order to obtain a quantitative understanding of this boost of correlation 
transmission by synchrony, in the following two subsections we will in turn investigate the mechanisms in the 
limit of low and high input correlations, respectively. 

2.3 Correlation Transmission in the Low Correlation Limit. For Gaussian white noise input and in 
the the limit of low input correlation, the correlation transmission can be understood by standard methods 
[44, 45, 35, 36]. This diffusion approximation assumes that the amplitudes of synaptic events are infinitesimally 
small. For uncorrelated Poisson processes and large number of afferents N, the theory is still a fairly good 
approximation. For small synaptic jumps approximate extensions are known [46, 47] and exact results can be 
obtained for jumps with exponentially distributed amplitudes [48]. However, in order to treat spiking synchrony 
in the common input to a pair of neurons, we need to extend the perturbative approach here. 

Before deriving an expression for the correlation transmission by a pair of neurons, we first describe the 
firing rate deflection of a neuron caused by a single additional synaptic impulse of amplitude J at t = on 
top of synaptic background noise. Within the diffusion approximation, the total afferent input to the neuron 
can be described by the first two moments p and a 2 (4). We determine the integral over the excursion of 
the firing rate h(t, J) = (s(t\ impulse of amplitude J at t = 0) — i> ut) caused by the additional impulse with 
respect to the base rate v out , as illustrated in Fig. 4B. An additional input with a stationary Poisson rate v 
has an effect r m Jv on the mean and ^r m J 2 v on the variance. The integral of the impulse response can be 

expressed as H(J) = f£°h(t,J) dt = dv °g} v) [47], where 

Voutty) ' s the equilibrium firing rate of the 

!/=0 

neuron receiving synaptic input causing a mean p + r m Jv and a variance a 2 + \r m J 2 v. Using the well known 
expression for the mean first passage time [49, 39] 

v-^{p,a) = T r + V^T m (F(y e )-F(y r )) (7) 
with 

F(y) = f V f(y)dy f(y) = e y 2 (etf(y) + l) 



and applying the chain rule we arrive at 



Ve - P _ V r - P 

V2a Vr ~ V2a 
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Figure 3: Correlation transmission of a pair of integrate-and-fire neurons. A Output spike synchrony as a 
function of input correlation p- in and for four different values of input synchrony p = 0, p = 0.001, p = 0.01 
and p = 0.1 (gray-coded). Dashed black line with slope 1 indicates p out = p ln . B Corresponding mean output 
firing rate of the neurons. C,D Cross-correlation functions at input correlations p m = 0.44 (C) and p m = 0.88 
(D) (indicated by dashed vertical lines in A) for the three values of input synchrony p as indicated in C. 



H(J) = / h(t, J) dt (8) 
Jo 

= aJ + (3 J 2 



with 



(^outr m ) 2 w|^ {f (y e ) - f(y r )) 



P = (foutTm) 2 ^^ {f(ye)ye - f{Vr)Vr) , 

where H is the expected number of additional spikes over baseline caused by the injected pulse. Assuming that 
the incoming events are pairwise correlated by a MIP [38] so that the sum of synchronously arriving events k 
is binomially distributed, we can expand (8) and set ,/ = kw, with k ~ B(cfN,p,k) to obtain 



H = wa^k- B(cfN,p,k) +w 2 fi^2k 2 ■ B(cfN,p,k) 

k k 

= waMi + w 2 /3M 2 , 



with Mi and Mi being the first and second moments of the binomial distribution 1 . 

If two neurons receive statistically the same input they fire with the same firing rate. Furthermore, due 
to the common input they both receive the same synchronous events with binomially distributed amplitudes 
J = kw. Any such event causes the same amount of excess spikes H(kw) in both neurons. If Si(t) denotes 
the spike train of neuron £=1,2, the cross covariance function is defined as 

K out (A) = ((si(A + t) - u ont )(s 2 (t) - iw))t, (9) 



1 The first four moments of the binomial distribution B(N,p) are 
Mi = Np, 

M 2 = Np(l - p + Np), 

M 3 = Np(l - 3p + 3Np + 2p 2 - 3JVp 2 + N 2 p 2 ) and 

M 4 = Np(l - 7p + 7Np + Yip 2 - 18Np 2 + 6N 2 p 2 - 6p 3 + UNp 3 - 6N 2 p 3 + N 3 p 3 ) 



8 



where the expectation value () t is taken over realizations of the input and over time t. K out (A) drops to zero 
for A — > oo. The only cause of covariation of the firing probability are the common events in the input. The 
sequence of events of amplitude J arriving with rate Aj contribute to the integrated covariance function as 



t (J) = J Kout (A,J)dA 

— oo 

oo T 

= / lim (2TXj)^- [ ( Sl (A + t)-v out )(s 2 (t)-v out )dtdA 

J T^oo 2T J 

-oo -T 

oo oo 

= Aj / (si(i) - v ont ) dt / (s 2 (t) - Vout) dt 



Aj / hit, J) dt I h{t, J) dt 





= \jH{J) 2 . 

The second step holds by exchanging the order of integrals and by inserting the expected number of events 
2TAj during time 2T. We now need to treat the excitatory and inhibitory input differently, since inhibitory 
events are not synchronized and furthermore additionally weighted by — g. Averaging over all common input 
events, the effective rate of an event of amplitude J = kw is given by \k w = ^ m B(cfN,p,k) we therefore 
obtain 

Sout = A m B (cfN,p, k) ■ H(kw) 2 + vc{\ - f)N ■ H(-gw) 2 , 

k 

and use H 2 (J) = a 2 J' 2 + 2a/3J 3 + (3 2 J 4 to obtain 

Kout = X m (a 2 w 2 M 2 + 2a/3w 3 M 3 + f3 2 w 4 lVh) (10) 
+Pc(l - f)N (a 2 (-gw) 2 + 2a^gw) 3 + P 2 (~gwf) , 

where M3 and M4 are the third and fourth moment of the binomial distribution 1 . 

In the limit of low input correlation we are interested in the output correlation on long time scales. In order 
to obtain this quantity, we need to normalize the integral of (10) by the integral of the auto-covariance of the 
neurons' spike trains. This integral equals FFj/ out [45], with the Fano factor FF. In the long time limit the 
Fano factor of a renewal process equals the squared coefficient of variation CV 2 [50], which can be calculated 
in the diffusion limit [39, App. Al]. Thus, we obtain 

Fig. 4A shows that the output spike correlation of a pair of neurons is fairly well approximated by p"^ t in 
the lower correlation regime. While the approximation is good over almost the whole displayed range of /?;„ 
for p = 0.001 and p = 0.01, for p = 0.1 the theory only works for values of pi n < 0.3 in agreement with 
previous studies [43, 36] applying a similar perturbative approach to the case of Gaussian input fluctuations. 

2.4 Correlation Transmission in the High Correlation Limit. In order to understand how the neurons are 
able to achieve a correlation coefficient larger than one, we need to take a closer look at the neural dynamics 
in the high correlation regime. We refer to the strong pulses caused by synchronous firing of numerous 
afferents as MIP events. Fig. 5A shows an example of the membrane potential time course that is driven by 
input in the high correlation regime. At sufficiently high synchrony as shown here, most MIP events elicit a 
spike in the neuron, whereas fluctuations due to the disjoint input alone are not able to drive the membrane 
potential above threshold. Thus, in between two MIP events the membrane potential distribution of each 
neuron evolves independently and fluctuations are caused by the disjoint input alone. Fig. 5B shows the time- 
dependent probability density of the membrane potential, triggered on the time of arrival of a MIP event. We 
observe that most MIP events cause an action potential, followed by the recharging of the membrane after it 
has been reset to V r at t = 0. After a short period of repolarization the membrane potential quickly reaches its 
steady state. The contribution of the cfN common, excitatory afferents to the membrane potential statistics 
is limited to those occasional strong depolarizations. Between two such events they neither contribute to the 
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Figure 4: Approximation of the output correlation in the limit of low input correlation. A Correlation 
transmission in the low input correlation limit. Data points show the output correlation Pout ms resulting from 
simulations, solid lines show the theoretical approximation p^ t (11). Dashed black line indicates p out = p- m . 
B Deflection of the firing rate with respect to base rate caused by an additional synaptic event at t = to- 



mean nor to the variance of V. Hence the effective mean and variance of the membrane potential are due 
to the disjoint input alone, given by p — Vq — g(l — f)cNv- m Fi and a 2 = [/(l — c) + g 2 (l — /)] Nv- m F2 
with F% = T m w and F2 — ^r m w 2 . Fig. 5C shows in gray the empirical distribution of the membrane 
potential between two MIP events after it has reached the steady state. It is well approximated by a Gaussian 
distribution with mean fl and variance a. The membrane potential can therefore be approximated as a 
threshold-free Ornstein-Uhlenbeck process [51, 52]. 

Let us now recapitulate these last thoughts in terms of a pair of neurons: In the regime of synchronized 
high input correlation (e.g. p = 0.1, p- m > 0.8), MIP events become strong enough so that most of them elicit 
a spike in both neurons. At the same time, the uncorrelated, disjoint sources (which can be considered as 
sources of noise) induce fluctuations of the membrane potential which are, however, not big enough to drive 
the membrane potential above threshold. Thus, while the input to both neurons still contains a considerable 
amount of independent noise, their output spike trains are (for sufficiently high pi n ) a perfect duplicate of the 
mother spike train that generates the MIP events in their common excitatory input, explaining the observed 
correlation transmission coefficient larger than unity. Note that this is the reason for the drastic decrease of 
the output firing rate in Fig. 3B, which in the limit of high input correlation approaches the adjusted input 
firing rate v- ln (Fig. 2B). 

We would like to obtain a qualitative assessment of the correlation transmission in the high correlation 
input regime. Since the probability of output spikes caused by the disjoint sources is vanishing, the firing 
due to MIP events inherits the Poisson statistics of the mother process. Consequently, the autocovariance 
function of each neurons' output spike train is a (5-function weighted by its rate v$ = A m Pi nst , where Pi ns t is 
the probability that a MIP event triggers an outgoing spike in one of the neurons. The output correlation can 
hence be approximated by the ratio 

Pout - (12) 

-*inst 

where P sync is the probability that a MIP event triggers an outgoing spike in both neurons at the same time. 
Note that the approximation (12) holds for arbitrary time scales, as the spike trains have Poisson statistics in 
this regime. In order to evaluate Pi ns t and P S ync, we use the simplifying assumption that the last MIP event 
at t = caused a reset of the neuron to V r — 0, so the distribution P(V,t) of the membrane potential evolves 
like an Ornstein-Uhlenbeck process as [52] 

P(V,t)= - 1 expf-^-f^H (13) 
with 

a 2 (t) =~a 2 (l-e-r 



which is the solution of (3) with initial condition V(0) = 0. We evaluate P; ns t from the probability mass of 
the voltage density shifted across threshold by an incoming MIP event as 
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Figure 5: Neural dynamics in the regime of high input correlation and strong synchrony. A Exemplary time 
course of a membrane potential driven by input containing strong, synchronous spike events. During the time 
period shown, five MIP events arrive (indicated by tick marks above Ve). The first four drive the membrane 
potential above the threshold Ve, after which V is reset to V r and the neuron emits a spike (dark gray tick 
marks above Ve). The fifth event is not able to deflect V above threshold (light gray) and the membrane 
potential quickly repolarizes towards its steady state mean ft (see text). B Time-resolved membrane potential 
probability density P(V, t) triggered on the occurrence of a MIP event at t = 0. Since most MIP events elicit a 
spike, after resetting V to V r the membrane potential quickly depolarizes and settles to a steady state Gaussian 
distribution. The slight shade of gray observable for small t just below the threshold Vg is caused by the small 
amount of MIP events that were not able to drive the membrane potential above threshold. C Probability 
density of the membrane potential in steady state. Theoretical approximation (black) was computed using p, 
and a (see text and (13)), empirical measurement (gray) was performed for t > 30 ms (gray dashed line in 
B). Simulation parameters were p = 0.1, c = 0.26 (c = 0.87) and v = 1.75 (v- m — 10) Hz. Other parameters 
as in Table 1. 



Pinst = J2 B ( 5 f N 'P> k ) ■ I dt X raS{t) ■ f dV P(V, t) , 

h — 1 



(14) 



Ve — kw 



where the survivor function S(t) = exp(— X m t) is the probability that after a MIP event occurred at t = 
the next one has not yet occurred at t > 0. So dtX m S(t) is the probability that no MIP event has occurred 
in [0, t] and it will occur in [t, t + dt] [50]. The binomial factor B is the probability for the amplitude of a 
MIP event to be kw and the last integral is the probability that a MIP event of amplitude kw causes an 
output spike [47]. We first express I(V,t) = JydV P(V,t) in terms of the error function using (13) with 



the substitution x — 



to obtain 



J(V,t) 



erf 



v e - m 

V^cr(t) 



where we used the definition of the error function erf(a;) 
integral in (14) with the substitution y 



e -A m t t0 



erf 



v-A(t) 



(15) 



2 



V2a(t) 

'= Jq e~ x2 dx. We further simplify the first 



o 

dy ( In y 



dt S(t) -I{V,t) = - -Z-yl[ V, 



Xy" V -A 



thus finally obtaining 



cfN 1 

Pinst = B ^f N i P^)- / HVe ~ kw, y) dy, (16) 

fc=i { 

where we introduced I(V,y) as a shorthand for (15) with jl(t) and a(t) expressed in terms of the substitution 
variable y as fi{y) = fx (l — y x mr m ^ anc | a(y) = a [\ — !/»"•'» J , following from (13). In order to approximate 
the probability P sync that the MIP event triggers a spike in both neurons we need to square the second integral in 
(14), because the voltages driven by disjoint input alone are independent, so their joint probability distribution 
factorizes, leading to 
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Figure 6: Approximation of the output synchrony in the limit of high input correlation. A Output spike 
synchrony as a function of input correlation in the limit of high input correlation and strong synchrony p = 0.1. 
Data points and solid lines show results from simulations and theoretical approximation (12), respectively. Gray 
code corresponds to the four different mean membrane potential values po as depicted in B, the input firing 
rate i/- m was 17.5 Hz, 13.7 Hz, 10.0 Hz and 7.2 Hz, correspondingly (from low to high po). The working point 
used in the previous sections corresponds to po = 10 mV, v- m = 10 Hz. The inset shows the output firing 
rate at the four working points. B Output spike synchrony as a function of the actual common input fraction 
c at the four working points. Dashed curves in A and B indicate p out = p m . 



Having defined Pi ns t and P syn c we now can calculate p out using (12) in the high input correlation regime. 
So far, we have considered both neurons operating at a fixed working point, defined by the mean and variance 
(4). Due to the non-linearity of the neurons we expect the effect of synchronous input events on their firing to 
depend on the choice of this working point. We therefore performed simulations and computed (2) using four 
different values for the mean membrane potential po (Fig. 6). This was achieved by an appropriate choice of 
a DC input current Io and accordingly adjusting the input firing rate V\ n in order to keep the mean firing rate 
constant (Fig. 6A, inset). The data points from simulations in Fig. 6A show that different working points of 
the neurons considerably alter the correlation transmission in the limit of high input correlation. At working 
points near the threshold (po = 11 mV) MIP events more easily lead to output spikes, thereby boosting 
the transmission of correlation, as compared to working points that are further away from the threshold 
(/io = 8 mV). Solid lines in Fig. 6A furthermore show that (12) indeed provides a good approximation of the 
output spike correlation when the input to both neurons is strongly synchronized. Obviously, the assumption 
has to hold that the probability density of the membrane potential is sufficiently far from the threshold, which 
for p — 0.1 is only the case if p ln > 0.75. Hence, the approximation becomes less accurate for lower input 
correlations, as expected. Note that, as opposed to Fig. IE, the effective common input fraction c in Fig. 6A 
is much lower than p m . Fig. 6B shows the same data as a function of the actual fraction of shared afferents 
c. It reveals that the gain of correlation transmission above unity is already reached at fractions of common 
input as low as c = 0.15 (for po = 11 mV), which is a physiologically plausible value. 

3 Discussion 

3.1 Summary of Results. In this work we investigate the correlation transmission by a neuron pair, using 
two different types of input spike correlations. One is caused solely by shared input - typically modeled as 
Gaussian white noise in previous studies [35, 36] - while in the other the spikes in the shared input may 
additionally arrive in synchrony. In order to shed light on the question whether cortical neurons operate as 
integrators or as coincidence detectors [18, 24, 25], we investigate their efficiency in detecting and transmitting 
spike correlations of either type. We showed that the presence of spike synchrony results in a substantial 
increase of correlation transmission, suggesting that synchrony is a prerequisite in explaining the experimentally 
observed excess spike synchrony [17, 21, 22], rather than being an epiphenomenon of firing rate due to common 
input given by convergent connectivity [8]. 

To model correlated spiking activity among the excitatory afferents in the input to a pair of neurons we 
employ the Multiple Interaction Process (MIP) [38], resulting in non-Gaussian fluctuations in the membrane 
potential of the receiving neurons. In this model the parameter p defines the pairwise correlation coefficient 
between each pair of N spike trains. If N is large enough and all spike trains are drawn independently (p = 0) 
the summation of all N spike trains is approximately equivalent to a Gaussian white noise process [40, 52]. 
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However, introducing spike correlations between the spike trains (p > 0) additionally allows for the modeling 
of non-Gaussian fluctuating inputs. Such correlations have a strong effect on the membrane potential statistics 
and the firing characteristics of the neurons. The fraction of common input c and the synchrony strength p 
each contribute to the total correlation between the inputs to both neurons. We show how to isolate and 
control the effect of input synchrony such that (1) a particular input correlation p ln can be realized by an 
(almost) arbitrary combination of input synchrony p and common input fraction c, and (2) the output firing 
rate of the neurons does not increase with p. This enables a fair comparison of transmission of correlation due 
to input synchrony and due to common input. We find that the non-linearity of the neuron model boosts the 
correlation transmission due to the strong fluctuations caused by the common source of synchronous events. 

Given a fixed input correlation, the correlation transmission increases with p. Most notably, this is the 
case although the output firing rate of the neurons does not increase and is for the most part constant, 
suggesting that the correlation susceptibility of neurons is not a function of rate alone, as previously suggested 
[35], but clearly depends on pairwise synchrony in the input ensemble. Previous studies have shown that the 
transmission of correlation of neuron pairs driven by Gaussian white noise can be approximated by employing 
Fokker-Planck theory and perturbation theory [44, 45, 35, 36]. In order to understand the effect of synchrony 
on the correlation transmission we extended this approximation to synaptic input of finite amplitudes. In the 
limit of low input correlation this extension indeed provides a good approximation of the output correlation 
caused by inputs containing spike correlations. Alternative models that provide analytical results are those of 
thresholded Gaussian models [53, 54] or random walk models [37]. In order to study transmission in networks 
with different architecture than the simple feed-forward models employed here, our results may be extended 
by techniques to study simple network motifs developed in [55]. 

Hitherto existing studies argue that neurons either loose correlation when they are in the fluctuation driven 
regime or at most are able to preserve the input correlation in the mean driven regime [56]. Here, we provide 
evidence for a qualitatively new mechanism which allows neurons to exhibit more output correlation than they 
receive in their input. Fig. 3A shows that in the regime of high input correlation the correlation transmission 
coefficient can exceed unity. This effect is observed at realistic values of pairwise correlations (p ~ 0.1) and 
common input fractions (c ~ 0.25). We provide a quantitative explanation of the mechanism that enables 
neurons to exhibit this behavior. We show that in this regime of high input correlation p- m the disjoint sources 
and the common inhibitory sources do not contribute to the firing of the neurons, but rather the neurons only 
fire due to the strong synchronous events in the common excitatory afferents. Based on this observation, in this 
contribution we derive an analytic approximation of the resulting output correlation beyond linear perturbation 
theory that is in good agreement with simulation results. 

3.2 Mechanism of Noise Suppression by Coincidence Detection. We presented a quantitative de- 
scription of the increased correlation transmission by synchronous input events for the leaky integrate-and-fire 
(LIF) model. Our analytical results explain earlier observations from a simulation study modeling synchrony 
by co-activation of a fixed fraction of the excitatory afferents [57]. However, the question remains what the 
essential features are that cause this effect. An even simpler model consisting of a pair of binary neurons 
is sufficient to qualitatively reproduce our findings and to demonstrate the generality of the phenomenon for 
non-linear units, allowing us to obtain a mechanistic understanding. In this model, whenever the summed 
input 7i 2 exceeds the threshold 6 the corresponding neuron is active (1) otherwise it is inactive (0). In Fig. 7 
we consider two different implementations of input correlation, one using solely Gaussian fluctuating common 
input (input G), the other representing afferent synchrony by a binary input common to both neurons (input 
S). The binary stochastic signal r)(t) has value A with probability q and otherwise, drawn independently for 
successive time bins. Background activity is modeled by independent Gaussian white noise in both scenarios. 
The input G corresponds to the simplified model presented in [35, cf. Fig. 4] that explains the dependence 
of the correlation transmission of the firing rate. In order to exclude this dependence, throughout Fig. 7 
we choose the parameters such that the mean activity of the neurons remains unchanged. As shown in the 
marginal distribution if the input current to a single neuron in Fig. 7B, in the scenario S the binary process 
r\ causes an additional peak with weight q centered around A. Equal activity in both scenarios requires a 
constant probability mass above threshold 9, which can be achieved by appropriate choice of as < <Jg- In 
scenario G the input correlation equals the fraction of shared input p in = c, as in [35], whereas in scenario S 
the input correlation is pi n = y^^+^r, where Var[r/] = q(l — q)A 2 is the variance of the binary input signal 
T)(t). Comparing both scenarios, in Fig. 7C-G we choose q such that the same input correlation is realized. 

As for our spiking model, Fig. 7C shows an increased correlation transmission due to input synchrony. 
This observation can be intuitively understood from the joint probability distribution of the inputs (Fig. 7D- 
G). Whenever any of the inputs exceeds the threshold [Ii t 2 > 0) the corresponding neuron becomes active, 
whenever both inputs exceed threshold at the same time (ii > 8 A 1% > 9), both neurons are synchronously 
active. Therefore, (/i) = J g dli J_ dli p(Ii, I2), the probability mass on the right side of 9 for input 
I\ (corresponding definition for (/z)), is a measure for the activity of the neurons. Analogously, (/1/2) = 
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Figure 7: Mechanistic model of enhanced correlation transmission by synchronous input events. A The 
detailed model discussed in the results section is simplified two-fold. 1) We consider binary neurons with a 
static non-linearity f(x) = H(x — 9). 2) We distinguish two representative scenarios with different models for 
the common input: G: Gaussian white noise with variance cxjq, representing the case without synchrony, or 
S: a binary stochastic process r/(t) with constant amplitude A, mimicking the synchronous arrival of synaptic 
events. In both scenarios in addition each neuron receives independent Gaussian input. B Marginal distribution 
of the total input I\$ to a single neuron for input G (gray) and S (black) and for p ln = 0.8. In input S 
the binary process rj alternates between (with probability 1 — q) and A (with probability q), resulting in a 
bimodal marginal distribution. The mean activity of one single neuron is given by the probability mass above 
threshold 8. We choose the variances Uq and er| of the disjoint Gaussian fluctuating input such that the mean 



Cov[/(Ji),/(J 2 )] 



as a function of the 



activity is the same in both scenarios. C Output correlation p out — 

^/Var[/(Ji)]Var[/(J 2 )] 

input correlation p- in (see A) between the total inputs 1\ and 1^. Probability q is chosen such that inputs G 
and S result in the same input correlation p ln . The four markers correspond to the panels D - G. D-G Joint 
probability density of the inputs I\, I2 to both neurons. For two different values of p; n the lower row (E,G) 
shows the scenario G, the upper row (D,F) the scenario S. Note that panel B is the projection of the joint 
densities in F and G to one axis. Brighter gray levels indicate higher probability density, same gray scale for 
all four panels. 



Ie° fe° dhdl2p(Ii, I2), the probability mass in the upper right quadrant above both thresholds is a measure for 
the output correlation between both neurons. Since by our model definition the mean activity of both neurons 
is kept constant, the masses and (/is) are equal in all four cases. However, the decisive difference between 
scenarios with inputs G and S is the proportion of (/1/2) on the total mass above threshold = (72) • This 
proportion is increased by the common synchronous events modeled as 77, observable by comparing Fig. 7D,E. 
The more this proportion approaches 1, (/1/2) — (fi) = {$2), the more the activity of both neurons is driven 
by 77 (Fig. 7F). At the same time the contribution of the disjoint fluctuations on the output activity is more and 
more suppressed. As the correlation coefficient relates the common to the total fluctuations, the correlation 
between the outputs can exceed the input correlation if the transmission of the common input becomes more 
reliable than the transmission of the disjoint input (cf. marker F in Fig. 7C). 

The situation illustrated in Fig. 7 is a caricature of signal transmission by a pair of neurons of a cell assembly. 
The signal of interest among the members of the assembly consists of synchronously arriving synaptic events 
from peer neurons of the same assembly. In our toy model such a volley is represented by an impulse of large 
amplitude A. The remaining inputs are functionally considered as noise and cause the dispersion of I\ and 
I2 observable in Fig. 7D-F. In the regime of sufficiently high synchrony (corresponding to large ^4) in Fig. 7F, 
the noise alone rarely causes the neurons to be activated, it is suppressed in the output signal due to the 
threshold. The synchrony coded signal, however, reliably activates both neurons, moving 1\ and I2 into the 
upper right quadrant. Thus a synchronous volley is always mapped to 1 in the output, irrespective of the 
fluctuations caused by the noise. In short, the non-linearity of neurons suppresses the noise in the input while 
reliably detecting and transmitting the signal. A similar effect of noise cancellation has recently been described 
to prolong memory life-time in chain-like feed forward structures [58]. 
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3.3 Limitations and Possible Extensions . Several aspects of this study need to be taken into account 
when relating the results to other studies and to biological systems. The multiple interaction process as a model 
for correlated neural activity might seem unrealistic at first sight. However, a similar correlation structure can 
easily be obtained from the activity of a population of N neurons. Imagine each of the neurons to receive 
a set of uncorrected afferents causing a certain mean membrane potential fj, and variance a 2 . The entire 
population is then described by a membrane potential distribution p(V). In addition, each neuron receives a 
synaptic input of amplitude w that is common to all neurons. Whenever this input carries a synaptic impulse, 
each of the N neurons in the population has a certain probability to emit a spike in direct response. The 
probability equals the amount of density shifted across threshold by the common synaptic event. Given the 
value p{Vg) and its slope §y(Vg) of the membrane potential density at threshold V$, the response probability 
is -fmst(w) = wp{Vg) — ^w 2 -^(Vg) + 0(w 3 ) to second order in the synaptic weight w. Employing the 
diffusion approximation to the leaky integrate-and-fire neuron, the density vanishes at threshold p(Vg) = 
and the slope is given by §&(Vg) — — ^ * 2" 1 [34]. The response probability hence is Pi ns t(w) — v °Xji m w 2 . 
For typical values of v out = 20 Hz, cr = 4 mV, and r m = 10 ms the estimate yields w = 5.7 mV to get the 
copy probability p = 0.1 used in the current study. Such a synaptic amplitude is well in the reported range for 
cortical connections [28]. As each of the neurons within the population responds independently, the resulting 
distribution of the elicited response spikes is binomial, as assumed by the MIP. Moreover, since our theory 
builds on top of the moments of the complexity distribution it can be extended to other processes introducing 
higher order spike correlations [57, 38]. 

The boost of output correlation by synchronous synaptic impulses relies on fast positive transients of the 
membrane potential and strong departures from the stationary state: An incoming packet of synaptic impulses 
brings the membrane potential over the threshold within short time. Qualitatively, we therefore expect similar 
results for short, but non-zero rise times of the synaptic currents. For long synaptic time constants compared 
to the neuronal dynamics, however, the instantaneous firing intensity follows the modulation of the synaptic 
current adiabatically [45, 59]. A similar increase of output synchrony in this case can only be achieved if the 
static f — I curve of the neuron has a significant convex non-linearity. 

The choice of the correlation measure is of importance when analyzing spike correlations. It has been 
pointed out recently that the time scale r on which spike correlations are measured is among the factors that 
can quantitatively change the results [3]. In particular, spike count correlations computed for time bins larger 
than the intrinsic time scale of spike synchrony can be an ambiguous estimate of input cross correlations [60]. 
Considering the exactly synchronous arrival of input events generated by the MIP, we chose to measure count 
correlations on a small time scale r = 1 ms (except for the analysis regarding the limit of low input correlation 
where we chose t = 100 ms). The observation that the difference between cross-correlation functions in 
absence and presence of input synchrony is localized to their peak values and is otherwise negligible reveals 
that this measure of zero-lag correlation fully captures synchrony-induced correlations. 

3.4 Conclusion. It had been proposed that the coordinated firing of cell assemblies provides a means for 
the binding of coherent stimulus features [14, 15, 16]. Member neurons of such functional assemblies are 
interpreted to encode the relevant information by synchronizing their spiking activity. Under this assumption 
the spike synchrony produced by the assembly can be considered as the signal and the remaining stochastic 
activity as background noise. In order for a downstream neuron to reliably convey and process the incoming 
signal received from the assembly, it is essential to detect the synchronous input events carrying the signal 
and to discern them from corrupting noise. Moreover, the processing of such a synchrony-based code must 
occur independently of the firing rate of the assembly members. We have shown that indeed the presence of 
afferent spike synchrony leads to increased correlation susceptibility compared to the transmission of shared 
input correlations. This finding, that the correlation susceptibility is not a function of the firing rate alone, 
demonstrates a limitation of the existing Gaussian white noise theory that fails to explain the qualitatively 
different correlation transmission due to synchrony [35]. We furthermore have shown that under realistic 
conditions cortical neurons are able to correlate their output stronger than the correlation they receive in 
their input, yielding a correlation transmission coefficient > 1. This observation is also accompanied by non- 
changing firing rates. This finding extends the prevailing view of neural correlations two-fold: The correlation 
susceptibility exceeding unity invalidates the description of correlation propagation as a 'transmission' per 
se. And moreover, the dependence of the output correlation on the type of the input correlation, not on its 
magnitude alone, demands an extended definition of a correlation transmission coefficient. We have shown 
in a mechanistic model how this coefficient exceeding unity results from the non-linearity of cortical neurons 
enabling them to actively suppress the noise in their input, thus sharpening the signal and improving the 
signal-to-noise ratio. This observation is in agreement with the mechanism of synfire activity where pulses of 
synchronized activity travel through feed forward structures in a stable manner [61]. From our findings we 
conclude that the boosting of correlation transmission renders input synchrony highly effective compared to 
shared input in causing closely time-locked output spikes in a task dependent and time modulated manner, as 
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observed in vivo [22]. 
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